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We study the characteristics of scattering processes at step edges on the surfaces of Strong Topo- 
logical Insulators (STI), arising from restrictions imposed on the S'-matrix solely by time reversal 
symmetry and translational invariance along the step edge. We show that the 'perfectly reflecting' 
step edge that may be defined with these restrictions allow modulations in the Local Density of 
States (LDOS) near the step edge to decay no slower than 1/x, where x is the distance from the 
step edge. This is faster than in 2D Electron Gases (2DEG) — where the LDOS decays as l/\/x 
— and shares the same cause as the suppression of backscattering in STI surface states. We also 
calculate the scattering at a delta function scattering potential and argue that generic step edges 
will produce a x~^^^ decay of LDOS oscillations. Experimental implications are also discussed. 

PACS numbers: 72.25.Dc, 73.20.-r, 75.30.Hx, 85.75.-d 



Strong topological Insulators (STI) are three- 
dimensional band insulators that have an odd number 
of gapless chiral modes on their surfaces [TH3]. These 
have recently been realized experimentally@]-[5] and are 
an active area of current research. The chiral states 
on the STI surfaces consist of time-reversed pairs of 
states propagating in opposite directions, between which 
backscattering by time-reversal invariant impurities and 
perturbations 5 is forbidden. Because of this, we can 
expect scattering at impurities or step edges on the 
STI surface to lead to outcomes that are substantially 
different from the case of the 2D Electron Gas (2DEG). 
This suppression of backscattering is also often quoted 
as a major reason for the topological protection of the 
STI surface states against surface impurities. While this 
suppression is precise for the case of exact backscatter- 
ing, it does not forbid scattering processes in which the 
particle gets reflected 'almost' backward - the formation 
of locahzed resonances/states near surface defects [71 |H] 
are still allowed. It is therefore instructive to test the 
response of these STI surface states to various surface 
defects and precisely quantify the restrictions imposed 
by the robust Time Reversal Symmetry (TRS). This 
knowledge will be an important ingredient in future 
attempts to design electronic devices based on these 
materials. Finally, the calculations in this paper also 
address a growing body of literature on the use of 
Scanning Tunneling Spectroscopy (STS) to investigate 
the role of isolated impurities and step edges on STI 
surfaces [3 IS [TU]. 

We shall address, in this paper, the scattering physics 
characteristic of a single set of chiral states - such as 
the Dirac states present on the surface of Bi2Te3 or 
Bi2Se3[31 E] — in the context of scattering from a step 
edge. We find the following: 

(i) TRS and unitarity impose a set of constraints on the 
reflection and transmission amplitudes (equation ([s])) ir- 
respective of the effective Hamiltonian describing the step 



LDOS near a perfectly reflecting step edge 




FIG. 1: (Color Online) Comparison between the LDOS near a 
'perfectly reflecting' step edge in the STI ( 16 1 and the 2DEG. 



Their envelopes, decaying as 1/x and respectively, are 

shown as dashed curves, k^i is the wave vector at energy uj. 



edge. We only require that the step edge does not violate 
TRS and is translationally invariant along its length, 
(ii) Suppressed backscattering leads to a substantial de- 
crease of LDOS modulation near the step edge on the 
surface of the STI — the LDOS is found to decay at 
least as fast as (Figure [T]) 
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In contrast, in 2DEGs this decay is slower ^ l/v^ [H]- 
(iii) We predict the existence of the 'perfectly reflecting' 
step edge (151 by using the scattering matrix restrictions. 



This perfectly reflecting wall produces LDOS modula- 
tions of the kind mentioned in ([T]) (shown in Figure [T]) . 
(iv) For a sharp step edge, which we approximate using 
a delta function, we have evaluated the reflection ampli- 



tude explicitly as a function of the potential strength ( 18 ) 
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and using this, find that for a 'strong' potential there are 
essentially no oscillations near the step edge. The LDOS 
decays monotonically as x~'^^'^ far from the edge, after 
an initial dip (Figure [S]) - similar to the long wavelength 
scattering observed in |10j. The x^'^^'^ decay law is also 
true for generic step edges due to TRS. 

We begin by formulating the general scattering matrix 
framework for the case of scattering at a surface step 
edge. We label the STI surface states by their band index 
s = ±1 and momentum fe = (kx-, ky)^ which together also 
determine their energies Eg^k- We choose i? = at the 
Dirac point and s = sgnEs^k- The defining characteristic 
of these states is that oppositely propagating states at 
a given energy are time-reversed partners (0 is the time 
reversal operator): 



|s, —k) cx Q |s, k) 



(2) 
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\r) = \s,k^,ky}R 



Step-edges don't usually break time-reversal symmetry 
unless, for e.g, they have embedded magnetic impurities. 
We can use this fact to find relations between the reflec- 
tion/transmission amplitudes defined above. For simplic- 
ity, we fix the phases of the incoming states by defining 
them to be the time-reversed versions of the appropriate 
outgoing states: 



-kx,ky)ji — Q\kx, ky) ji 
\kxTky)L ^ 0| kx, ky) 



(4) 



Using this convention and the relation 0^ = — 1, we can 
prove the following constraints that are similar to the 
Stokes' relations in Optics: 

vikx) ky) — f(^kx, ky), r {kx, ky) — r {kx^ ky") 

t{kx,ky) = t {kx, ky) (5) 

As an example, we prove the first relation explicitly: 




kx^ ky) 



R 



skx > 



FIG. 2: (Color Online) The scattering problem at a step edge, 
at an energy Es,k and a j/-momentum ky. The green (incom- 
ing) and dashed violet (outgoing) arrows label the 'incoming' 
and 'outgoing' states respectively. 



We assume that there is translational invariance in the 
direction parallel to the step edge, say y (See Figure |2]). 
Thus, the step edge can be represented by a finite-range 
modification to the Hamiltonian that preserves the y- 
moment um 1 1 11 [T^ . In some cases, it may be modeled as 
an electrostatic potential V{x) that is a function of only 
the perpendicular coordinate x. Scattering processes at 
the step leave ky - the y-momentum ~ unchanged. In 
what follows, we shall first embody the effect of the step 
edge in a scattering matrix (the S'-matrix) . We label the 
scattering process as in Figure [2] and define the reflec- 
tion and transmission coefficients to be (r, r') and (i,t') 
respectively. The energy of the scattering process is ■ 
The states are taken to belong to the s"^ band - this 
is not mentioned explicitly, for brevity, in the following 
equations. RjL denotes the half-space (right/left) these 
states belong to. The 5-matrix is 
vijzx^ky) t (^kxjky 



ky) — 



t[kx,ky) r [kx,ky) 



(3) 



B.{kxTky\S\ kx,ky)j^ B.{kx,ky\S\kx,ky) ij 
l{ kx, ky\S\ kx, ky) n l{ kx , ky\S\kx, ky) j 



r(kx,ky) — [kx,ky\S\ kx,ky)ji 
=R {Q{-kx,~ky)\S\ - kx,ky)R 

{e{-kx,ky)\S\Q^(-kx,~ky))R 

= R{kx, kx, ky) R = —r{kx, ky) 



(6) 



The antiunitary nature of the time reversal operator was 
used in the third line and = —1, corresponding to 
spin-1/2 states, was used in the last step. 

The unitary nature of the scattering process requires 
SS'f = I and so: 



|rp + |t|2^|/|2 + |t'|2 = i 
r*t' + r't* = 
(^ |r|2 + |t'|2 ^ |/|2 + |i|2 ^ 1 also) 



(7) 



If we are able to specify S, the asymptotic forms of the 
new energy eigenstates, i.e, the wavefunctions outside the 
influence of V{x), are given by: 



\k)T = 



kx, ky) R + r(fc) \kx, ky) + t{k) \ kx, ky) j 



V2 



i.Nncw _ \kx,ky)L + r'{k)\ kx,ky)^+t'{k)\kx,ky)j^ 
I'^/L — 



V2 



(8) 



The normalizations follow from ([7|. Since we have these 
new eigenstates of the Hamiltonian, it is straightforward 
to calculate the modified LDOS away from the step edge: 



3 



p{ijj, cc > 0) = 



TT 



27r 



s,k) 



dkxdky 
^Re 



(9) 



Here, p*^"^ is the LDOS in the absence of the step and 
simphfications have been made using time reversal sym- 
metry. At this stage we can approximate the band struc- 
ture near the Dirac point as that arising from the Dirac 
hamiltonian 1161 



y. — vcr ■ k 



(10) 



Here, <t is the actual spin operator, up to a rotation. 
The outgoing eigenstates of this Hamiltonian are given 
by {skx ^ for the R/L cases): 
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ih-r 



, Es,k = svk (11) 
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The phases of the incoming eigenstates are set by Q 

( 
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We have used 6 = iayK, K being the complex conjuga- 
tion operator, in the above equation. Using this setup, 
([9]) can be simplified to yield (for a; > 0): 



p«')(w) 



2sgn^ Tj'^'gi^^Ij^ ^^^(^)g2^fc..cost 



(13) 



Here, p'"-'(cli) = |cLi|/(27n;^), k^^ — lj/v is the momen- 
tum at energy uj and the reflection amplitude r^{6k) = 
r(sgna;, fe^^) needs to obey the property ^ 



(14) 



We now note that the 'perfect' reflector needs to have 
the form (no transmittance except at normal incidence) 



(0) 



isgn(0) 



(15) 



which tells us that the LDOS becomes (a; > 0): 



5p{u},x) 2sgnaj 



tt/2 



{cos{2ki^x + Suj) — cos S, 
agnu irki^x 



dO sin 6* sin {2k^x + 6^^^\e\) 



(when 5 is 0-indep.) 



(16) 



In Figure[TJ these LDOS modulations are compared with 
those near a perfectly reflecting step edge in a 2DEG. We 



see that even for 'perfect' reflection the step edge is only 
able to create LDOS oscillations that decay as l/x, as 
opposed to the slower 1/y/x decay in the 2DEG casejllj. 
In a recent paper (Figures 3a-e in [ICF) oscillations of 



the form given by ( 16 1 have been observed at energies far 
from the Dirac point where the wavefunction wavelength 
is comparable to or shorter than the width of the step 
edge. However, in that region the band surface exhibits 
hexagonal warping |13| and the dominant scattering sig- 
nature comes from the scatterings between the states at 
the adjacent Af-points (the step edge is oriented per- 
pendicular to the TM direction). These do not involve 
the special processes near perfect backscattering in the 
STI that we have considered above. The 1/x behavior 
arises simply because the bands are not perpendicular to 
the TK direction near the hexagon vertices[17]. In Fig- 

LDOS near a delta-function step edge 
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Case U = -3 
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Case U = -25 
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FIG. 3: (Color Online) LDOS above the Dirac point near a 
strong attractive (5-function step edge potential (171 for vari- 



ous values of U (solid curves with different colors and thick- 
nesses). Notice how the oscillations get damped for large U . 
The dashed curve is the envelope for the 'perfectly reflecting' 
step edge on the STI (Figure [T]), shown here for comparison. 



ures 3f-h of the above-mentioned paper [TD] the authors 
also observed that at lower energies, closer to the Dirac 
point, the character of these LDOS modulations change 
qualitatively. They are much less pronounced exhibit an 
almost oscillation-less decay. We propose to explain this 
behavior by modeling the step edge as a delta-function 
potential 



V{x) = U6{x) 



(17) 
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when considering the scattering of long wavelength states 
near the Dirac point. Using this with the original Hamil- 



tonian (10), we obtain the reflection amplitude|14j: 



AU sm{e) 



[U^ ~ 4) cos(6l) - 4isgnwf7 



(18) 



Note that the reflection amplitude is antisymmetric 
about 9 = 0, and is also linear around that region. The 
LDOS obtained for this delta- function potential is shown 
in Figure [3] and displays minimal oscillations when U is 
large and negative and E > Q. There is a dip in the 
LDOS near the wall followed by an almost oscillation- 
less decay (for large potential strength) away from the 
wall - similar to what was observed in [TU]. The dip 
can be explained by the existence of bound states near 
the wall: for U < —2 these occur at negative energies 
but the transferred spectral weight leads to the dip at 
positive energies. 

The long-distance decay of the oscillations is given by 



-3/2 



which arises out of the proportionality between 



r{9) and 9 near 9 = 0. Because of the antisymmetry of 



r{9) about 6* = ( 14), for generic step edges we will also 
have r{9) oc 9 near 9 = and this will result in the almost 
coherent scatterings near 9 = interfering far away from 
the edge to yield a spatial decay law of a;^^/^[15 . This 
power law may be obtained by considering the scaling 
behavior of the integrand in (13), near 9 = 0. 



It is possible that for high energies, near a step edge 
perpendicular to the TM direction in Bi2Te3 (or with any 
orientation in Bi2Se3 which has no hexagonal warping), 
one will observe a 1/x decay of oscillations at the wave- 
vector equal to the diameter of the band surface along 
the K — T — K direction if the step edge is perfectly non- 
reflecting - a quality that possibly arises due to roughness 
along its normal direction. Generically, due to the scaling 
relation r{9) oc 9 near 9 = 0, one should be able to 
observe a faster x^"^/^ decay of the oscillations. 

In summary, we have used the chiral nature of the STI 
surface states and time reversal invariance to impose re- 
strictions on the reflection and transmission amplitudes 
of the scattering process at a step edge on the surface of 
a Strong Topological Insulator. This allowed us to define 
the 'perfect reflecting step edge' for these chiral surface 
states and we found that the amplitude of LDOS rip- 
ples caused by such an edge decay no slower than 1/x — 
faster than in the analogous case for a 2DEG, where the 
decay occurs as For a generic smooth step edge, 

the decay law is x^^^'^. We give possible reasons for the 
LDOS features seen in experiments [TU] — damped oscil- 
lations and a monotonic decay of the LDOS modulation 
far from the step edge for small energies and a l/x decay 
at higher energies - by arguing that the scattering of long 
wavelength states near the Dirac point occurs due to an 



the band surface yields the l/x decay at higher energies. 
To the best of our knowledge, this is the flrst application 
of this general ^'-matrix formalism (with the imposition 
of TRS) to the scattering of STI surface states. 
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